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ABSTRACT 

We report the discovery of a star-forming loop around the young, Crab-hke supernova remnant 
(SNR) G54. 1-1-0.3 using the AKARI infrared sateUite. The loop consists of at least 11 young stellar 
objects (YSOs) embedded in a ring-like diffuse emission of radius ^ 1'. The YSOs are bright in the 
mid-infrared and are also visible in the Spitzer Space Telescope Galactic plane survey images. Their 
Spitzer colors are similar to those of Class II YSOs in [3.6] — [5.8] but significantly redder in [8] — [24], 
i.e., <[3.6] — [5.8]< 1.2 and 5 <[8] — [24]< 9. Most of them have near- infrared counterparts in the 
2MASS JHKs images, and some of them have an optical counterpart too. Their JHKg colors and 
magnitudes indicate that the YSOs are massive ( ^ 10 Mq) pre-main-sequence stars at the same 
distance to the SNR, i.e., 8 kpc, which supports the association of the star-forming loop with the 
SNR. The dereddened spectral energy distributions are similar to early Herbig Be stars, which are 
early B type pre-main-sequence stars with inner disks that have been destroyed. The confinement to 
a loop structure indicates that the YSOs are young, i.e., ^ 2 Myr. We propose that their formation 
is triggered by the progenitor star of G54. 1-1-0.3, which has a mass of ^ 15 AIq. The triggering must 
have occurred near the end of the progenitor's life, possibly after it had evolved off the main sequence. 
Subject headings: ISM: individual (G54. 1-1-0.3) — infrared: stars — stars: formation — stars: pre- 
main-sequence — supernovae: general — supernova remnants 



1. INTRODUCTION 

Massive stars can trigger the formation of a next gen- 
eration of stars. Their strong stellar ultraviolet radia- 
tion and winds can compress the ambient medium so 
that it becomes gravitationally unstable and collapses. 
Numerous examples supporting the scenario of triggered 
star formation have been proposed, from young stars 
formed in small bright-rimmed globules inside HII re- 
gion s to OB as sociations around kiloparsec supershells 
(see lElmegreen 2002; Zinncckcr & Yorkc 2007, and ref- 
erences therein). 

In this Letter, we report the discovery of a unique 
example of triggered star formation, namely trigger- 
ing by the progenitor star of a young, core-collapse 
SNR. The SNR (G54.H-0.3) is Crab-hke, with centrally- 
brightened synchrotron e r nission in radi o and X-rays 
(|Velusamv fc Beckeil Il988t iLu et al.l |2002[ ) . The extent 
of the remnant is 120" x 80" and is much fainter than 
the Crab, e.g., 0.5 Jy versus 1040 Jy at 1 GHz. At the 
center, a 136 ms radio/X-ray pulsar with a charact eristic 
age of 2900 yr has been discovered (|Lu et al.ll20"ol ). The 

^ Department of Physics and Astronomy, Seoul National 
University, Seoul 151-742, Korea; koo@astrohi.snu.ac.kr, 
hglee@astro.snu.ac.kr, sshong@astro.snu.ac.kr 

^ Department of Physics and Astronomy, University of Califor- 
nia, Berkeley, CA94720, USA; cmckee@astro.berkeley.edu 

^ Astronomy & Astrophysics Department, Pennsylvania State 
University, University Park PA16802, USA; lee@psu.edu 

Department of Astronomy and Space Science, Sejong Univer- 
sity, Seoul 143-747, Korea; jelee@sejong.ac.kr 

^ Department of Astronomy & Astrophysics, University of 
Toronto, Toronto ON M5S 3H4, Canada; moon@astro.utoronto.ca 
Institute of Space and Astronautical Science, JAXA, Kana- 
gawa 229-8510, Japan; kaneda@ir.isas.jaxa.jp 

''' Department of Astronomy, University of Tokyo, Bunkyo-ku, 
Tokyo 113-0003, Japan; onaka@astron.s.u-tokyo. ac.jp 



distance to the SNR is 8 kpc (§ 4). This object provides 
a laboratory in which to explore the physics of triggered 
star formation. 

2. DISCOVERY OF THE IR LOOP 

A mid-infrared (MIR) observation of G54. 1-1-0.3 was 
done using the Infrared Camera (IRC) aboard AKARI 
on 2007 April 17. The observation was performed in a 
framework of the AKARI mission program: the ISM 
in our Galaxy and nearby galaxies (ID: 1401070.1, PI: 
H.K.) The IRC is equipped with three wave band chan- 
nels covering the 2.6-26 fim wavelength range w ith a 
10' X 10' field-of-view (FOV) (lOnaka et al.l [2001 . We 
used the MIR-L channel in IRC02 mode, which gave two 
band images, L15 and L24, centered at wavelengths of 
15.58 and 22.89 /im, respectively. Their angular resolu- 
tions are 5. "7 (L15) and 6. "8 (L24). The total on-source 
integration time was 196 s for each band. The basic 
calibration and data handling were processed using the 
standard IRC Imaging Data Reduction Pipeline version 
070104 ^. The positional uncertainty is 0."15 at the 1-cr 
confidence level. 

Fig. 1 shows the AKARI L15 image of G54.1-h0.3 
with the 4.85 GHz radio continuum image overlaid. The 
radio image was made from archive data of the Very 
Large Array (VLA) of the National Radio Astronomy 
Observatory ^; the data were originally obtained by 
iVelusamv fc Becked (fT98l . The AKARI image shows 
a well-defined, bright, rin g-like struct ure at the posi- 
tion of the SNR (see also isiand l2007l ). The ring ap- 



° http:/ /www. ir. isas.jaxa.jp/ASTRO- 

F/Observation/DataReduction/IRC/ 

^ The National Radio Astronomy Observatory is a facility of 
the National Science Foundation operated under cooperative agree- 
ment by Associated Universities, Inc. 
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pears partially complete with its northeastern portion 
opened, with faint ridges extending out from the ends 
of the bright ring. The bright portion of the ring (here- 
after the "main ring" ) is elongated along the northwest- 
southeast direction and has an extent of ~ 105" x 54", or 
4.1 X 2.1 pc at 8 kpc. The ring is composed of both diffuse 
and compact sources. Most of the compact sources are 
point-like except the bright one at the northwestern end. 
There are more than 10 point sources distributed along 
the main ring and several others at the 'tips' of the faint 
ridges. Most of them are bright in the 24 /im image and 
have colors quite different from foreground/background 
stars. This can be seen in Fig. 1 (right), which is a 
three-color image produced from the Spitzer IRAC 5.8 
urn (B), AKARI 15 /^m (G), and Spitzer MIPS 24 /im 
(R) images. We obtained the Spitzer image of the source 
from the GLIMPSE (Galactic Legacy Infrared Mid-Plane 
Survey Extraordinaire, Release II) and MIPSGAL legacy 
projects. There are 11 sources with colors significantly 
redder than the other stars in the field (see also the inset 
in Fig. 2). 

3. YOUNG STELLAR OBJECTS IN THE IR LOOP 

We conducted photometry of a 4.'5-radius area sur- 
rounding the IR loop using the AKARI 15 and 24 fim 
and the Spitzer MIPS 24 /im and 70 /im images. We 
have identified 151 sources detected in at least one of 
the 15 or 24 /im bands; most of them were detected at 
15 /tm while 76% were detected at 24 /im. At 70 /tm, 
only one source was identifiable and only upper limits 
are derived for the rest. The photometry of stars in the 
IR loop is not straightforward because of the bright dif- 
fuse emission from the loop. We did PSF photometry to 
minimize the contamination due to the diffuse emission. 
Sources are crossidentified with the Spitzer GLIMPSE 
catalog, which contains fluxes of the four IRAC bands, 
and also with the 2MASS Second Incremental Release 
Point Source Catalog. 

Fig. 2 shows the Spitzer [3.6] -[5.8] versus [8] -[24] 
color-color plot of the 64 sources detected in the IRAC 
and MIPS 24 /im bands with a [24] magnitude error of 
less than 0.2 mag. The stars clustered near (0.0, 0.0) are 
probably main-sequence or giant stars along the line of 
sight. The sources along the IR loop are well defined in 
this diagram by their strong [8] — [24] excess and moder- 
ate or small [3.6] — [5.8] excess, i.e., 5 <[8] — [24]< 9 and 
<[3.6] — [5.8]< 1.2. The sources within this box are 
marked by filled circles. The half-filled circles with ar- 
rows in the box represent the sources not detected in one 
of the IRAC 5.8 and 8.0 /im bands but are included be- 
cause they are likely to fall into the box based on their 
MIR colors. We marked their positions using a shorter 
IRAC band, e.g., if a source is not detected at 8.0 /im 
we used [5.8] instead of [8.0] and placed an arrow. The 
inset shows that all sources in the box are distributed 
along the IR loop, although the loop of YSOs is quite dis- 
torted. Their colors are different from those of previously 
studied YSOs in nearby star-forming regions. We show 
this by drawing boxes (dotted lines) in Fig. 2 represent- 
ing the areas occupied by young embedded protostars 
with both circumstellar disks and envelopes (Class O/I) 
and pre-main sequence stars with significant circumstel- 
lar disks (Class II) de termined by previous studies (e.g., 
iMuzerolle et al.ll2004D . It is clear that the YSOs in the 



IR loop do not fall into the usual Class O/I or Class II 
areas and that this cannot be due to extinction. There 
are many sources that fall into and around Class II areas. 
They are outside the IR loop and will not be discussed 
further in this paper. 

Fig. 3 (bottom) shows a J — H versus H — Kg color- 
color diagram for all 110 MIR sources identified in the 
2MASS JHKg bands. The small dots represent sources 
not detected in Spzteer MIPS 24 /im or IRAC bands while 
the other symbols have the same meaning as in Fig. 
2. There are only nine filled/half- filled circles because 
two YSOs with large [3.6] -[5.8] (Fig. 2) are not iden- 
tified in the 2MASS bands. The dashed line represents 
the m ain-sequence (MS), which is from iBessel fc Brettj 
(|1988f ) for stars later than AO and is based on calculations 
of the JHKg colors fro m Kurucz model atmospheres 
( Castelh fc Kurucj[2003l ) for stars earher than AO. The 
dotted line represents interstellar redde ning line of an 
AO st ar following the reddening law of iBessel fc Brett! 
(|1988f l. Note that the YSOs are clustered around the 
positions of reddened OB stars. Another way to see this 
is from the color-magnitude {Kg vs. H — Kg) diagram 
(Fig. 3 top). In this diagram, the dashed line represents 
the MS at 8 kpc calculated from the Kurucz model atmo- 
spheres and using the relation between effective temper- 
ature a nd radius of MS star s. We adopt the relation of 
iMartin s. Schaercr, & Hillicr (l2005l ) for 0-type stars and 
that of I Schmidt-Kaler (1982i) for stars later than BO. For 
BO stars, we interpolate the two with an adopted tem- 
perature of Toff = 30, 000 K. The dotted lines represent 
the reddening lines that bound the YSOs. The YSOs are 
again clustered around the positions of reddened early- 
type stars. The YSOs are not aligned along a single MS 
line, but that may be due to slightly different amounts 
of extinction toward individual sources. The extinction 
toward individual objects can be derived if we assume 
that they are on the MS, and it ranges from Ay = 6.9 
to 9.2 with a mean of Ay = 8.0 ± 0.7 mag. The mean 
extinction agrees very well with the extinction implied 
from the X-ray absorbing H columns to the SNR, i.e.. 
Ay = 8.4 ±0.5 mag from N{H) ^ 1.6(±0.1) x 10^^ cm'^ 
(|Lu et al.l l2002) . We conclude that the JHKg color-color 
and color-magnitude diagrams are consistent with the 
YSOs being massive pre-main-sequence stars at 8 kpc. 
The spectral types determined from their Kg magnitudes 
in Fig. 3 (top) are B1.5 to 08, which have Toff =23,000 
- 35,000 K, masses of 12-21 A/0, and luminosities of 
(1-10) xlO^ Lq. 

Fig. 4 shows the SED of the most luminous source 
in the Kg band (the brightest one in the southern part 
of the loop in Fig. 1) as an example of the YSO 
SEDs. The filled symbols represent dereddened fluxes 
using the value Ay = 6.9 mag derived from Fig. 3 
and a n interstellar reddening law with Ry = 3.1 (jPraind 
I2OO3I ). The dereddened SED is characterized by a star- 
like SED with a MIR/FIR excess at 8 /im, as we 
could have inferred from Figs. 2 and 3. The spec- 
trum at ^ 5 /im matches reasonably well the SED of 
a 35,000 K MS star, the temperature of which was deter- 
mined from its Kg magnitude (Fig. 3). This spectrum 
is quite similar to the SED of early Herbig Be (HBe) 
stars (spectral type B0-B5; group III ) ("Hillenb rand et al] 
ll992tlMalfait. Bogaert. fc Waelkensiri998 1. A dip at 6- 
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10 /im similar to that in Fig. 4 was observed in HBe 
stars and has been interpreted as a physical hole in the 
dust distribution caused by the brea k-up of inner disks 
(|Malfait. Boeaert. fc Waelke^lTMsh . If we consider the 
MIR/FIR excess at ^ 8 fim to arise from an optically 
thick disk with a constant effective temperature and fit 
the SED with a 35,000 K star plus a blackbody, we ob- 
tain a disk with a temperature of 190 K and a lumi- 
nosity of 160 Lq, although this fit is not quite satis- 
factory. Most YSOs in the IR loop have SEDs similar 
to Fig. 4, with some objects showing weak MIR excess 
also. The ratio of the disk to stellar luminosity is ^ 
1%, showing that there is not much circumstellar mate- 
ri al near the YSOs . If w e use the empirical relationship 
of lNakai ima et all (|2005[ ) between stellar masses and K 
magnitudes for group III HBe stars, the masses of the 
YSOs range from 12 to 18 Mq which agrees with our 
estimates above. The 2MASS survey has a limiting Kg- 
band sensitivity of 14.3 mag, corresponding to a 6 AIq 
star at 8 kpc with 8 mag of extinction; no additional 
MIR sources are seen between 13.2 and 14.3 mag in the 
vicinity of the IR loop, however, which suggests a non- 
standard IMF for this region. 

4. SNR/IR LOOP ASSOCIATION AND THEIR DISTANCE 

The IR loop surrounds the X-ray /radio bright central 
portion of the SNR almost perfectly. The SNR is some- 
what extended along the northeast-southwest direction, 
suggesting that its expansion has been blocked toward 
the IR ring. We note that the radio contours are de- 
formed around the northwest, which is consistent with 
encountering a dense medium there. The positional co- 
incidence with a detailed spatial correlation strongly sug- 
gests that the IR loop and the SNR O54.H-0.3 are as- 
sociated. The agreement between the extinction toward 
the IR loop and the X-ray absorbing column toward the 
SNR further supports the association (§ 3). The absence 
of direct evidence for the interaction of the SNR shock 
with the dense gas may be understood if we assume that 
the IR loop is a partial shell in a low density medium, 
so that the SNR shock has been able to propagate well 
beyond the loop now. 

The distan c e to 54. 1-1-0. 3 has been determined by 
lOamilo et al.l ()2002f ). who derived the dispersion mea- 
sure of the central pulsar (308 cm~^ pc) and noted that 
it corresponded to a distance of 8-12 kpc, depending 
on the electron-distribution model of the Oalaxy. We 
have obtained an Hi absorption spectrum of the SNR 
using the VLA Oalactic Plane Survey data (|Stil et al.l 
l2006l i. The Hi spectrum of the SNR shows Hi ab- 
sorption over all positive velocities, which indicates that 
054. 1-1-0.3 is beyond the tangent point (> 5 kpc). The 
absorption peak occurs at -1-23 km s^^, and at this ve- 
locity there is a faint ^^00 J = 1-0 emission coinci- 
dent with the IR loop in the Boston-University-Five Col- 
lege Radio Astronom y Observatory Oalactic Ring Survey 
(| Jackson et al.ll2006f ). The emission is weak (T^ < 0.7 K) 
and implies a mass ^ 400 Mq. The kinematic distance 
correspo nding to -1-23 km s~^ i s 8.2 kpc using the rotation 
curve of iBrand fc Blitj (|1993[ ). who adopt Rq = 8.5 kpc 
and Oq = 220 km s~^. This is consistent with the dis- 
tance derived from the pulsar DM. We adopt 8 kpc as 
the distance to 054. 1+0.3 and the IR loop, and note 
that this is the distance to the Perseus spiral arm in this 



direction. 

5. POST MAIN-SEQUENCE TRIGGERED STAR 
FORMATION? 

The alignment of the YSOs along a small, loop-like 
structure strongly suggests that we are observing star 
formation triggered by some mechanism originating at 
the center. It cannot be the SNR 054. 1+0.3, however, 
because the remnant is only a few thousand years old. 
Instead, the location of the 054. 1+0.3 pulsar near the 
center of the IR loop indicates that it is most likely the 
progenitor star of the SNR and that the massive pro- 
genitor star produced an HII region/ wind bubble that 
compressed the ambient medium into a shell that frag- 
mented and collapsed to form the YSOs. 

What is more interesting about this system is that the 
triggered star formation must have occurred near the end 
of the progenitor's life, possibly when it was in its post 
main-sequence phase. The ages of early Herbig Be stars 
are not well known because they usually cannot be esti- 
mated from the H-R diagram (e.g. , ,Testi, Palla, fc Natta| 
Il998f ). The surrounding diffuse emission in the IR loop, 
which might be the emission from dust grains in natal 
cloud heated either by the YSOs or/and other nearby B 
stars, indicates that they are young. We can set a limit 
on the age of the YSOs by noting that 100 Mq of stars 
in a region with a radius of 2 pc should have a ID veloc- 
ity dispersion of about 0.2 km s~^. In order for these 
stars to have a regular spacing of about 0.7 pc as ob- 
served, they could not have moved more than ^ 0.35 pc, 
which sets a limit on their age of about 2 Myr. This 
limit is significantly shorter t han the main-seq uence life- 
time of the progenitor star: IChevalierl ()2005[ ) identified 
O54.1+0.3 as the remnant of Type IIP SN, with a pro- 
genitor mass ~ 10 — 25M0. The destructive impact of 
a massive star on the amb ient molecular c loud increases 
rapidly with stellar mass (jOhevalieil 1 1999t ) . so the small 
size of the IR loop and the regular structure seen in the 
spatial distribution of the YSOs suggest that the progen- 
itor was near the lower end of this mass range. A 15M0 
star has a lifetime ~ 13 Myr, and it spend s about 10% of 
its lif e in the post main-sequence phase (jSchaller et al.l 
119921 ). The post main-sequence lifetime increases for 
lower stellar masses. Therefore, it is possible that the 
YSOs formed during the post main-sequence phase, al- 
though it is likely that the progenitor had a significant 
dynamical effect on the natal cloud while it was still on 
the main sequence. 

We do not have enough information to determine 
why the triggering was delayed, i.e., why it did not 
occur durin g the first few Myr as theories predict 
(e.g., iHosokawa fc Inutsukal l2006( l. Two possibilities 
can be ruled out: The triggering was most likely not 
due to the wind associated with the red supcrgiant 
(RSO) phase of evolution, since a typical RSO wind 
is far too weak to have a sign ificant effect on the shell 
()Chevalier fc Emmerin^ I1989D . In any case, the RSO 
wind in 054. 1+0.3 may have been particularly weak 
since no SNR shock propagating into ambient or cir - 
cumstellar medium has been detected () Ohevalieij 1 2 5() . 
One might also think that the temperature and ioniza- 
tion of the ambient molecular gas would drop signif- 
icantly whe n the progeni t or ev olved off the main se- 
quence, but ISpaans et al.l ()1994[ ) showed that the tem- 
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perature and ionization in molecular gas that is near a 
star do not vary much for effective temperatures in the 
range (6 — 30) x 10'^ K. A potential explanation for the 
delay in the triggering of the star formation is that the 
natal cloud was shielded by other material until a few 
Myr ago, when it was exposed to the star by a combi- 
nation of motion of the star relative to the cloud and 
photoevaporation of the intervening gas. Future obser- 
vations will reveal the properties of the newly formed 
stars and shed light on the triggering of star formation 
by evolved stars. 

This work is based in part on observations with AKARI 



, a JAXA project with the participation of ESA, and also 
in part on observations made with the Spitzer Space Tele- 
scope, which is operated by the Jet Propulsion Labora- 
tory, California Institute of Technology under a contract 
with NASA. We wish to thank David Hollenbach and 
Stan Kurtz for their helpful comments. This work was 
supported by the Korea Science and Engineering Foun- 
dation (ROl-2007-000-20336-0) and the Korea Research 
Foundation (R14-2002-058-01003-0). The research of 
CFM is supported by the National Science Foundation 
through grants AST-0606831 and PHY05-51164. This 
publication makes use of molecular line data from the 
Boston University-FCRAO Galactic Ring Survey (GRS). 



REFERENCES 



Bessell, M. S. & Brett, J. M. 1988, PASP, 100,1134 
Brand, J., & Blitz, L. 1993, A&A, 275, 67 

Camilo, F., Lorimer, D. R., Bhat, N. D. R., Gotthelf, E. V., 
Halpcrn, J. P., Wang, Q. D., Lu, F. J., & Mirabal, N. 2002, 
ApJ, L71 

Castelli, F. & Kurucz, R. L. 2003, in lAU Symp. Ser. 210 edited 

by N. Piskunov, W.W. Weiss, & D.F. Gray, 20 
Chevalier, R. A. 1999, ApJ, 511, 798 
Chevalier, R. A. 2005, ApJ, 619, 839 

ChevaUer, R. A., & Emmering, R. T. 1989, ApJ, 342, L75 
Draine, B. T. 2003, ARAA, 41, 241 

Elmegreen, B. G. 2002a, in lAU Symp. Ser. 207 edited by D. 

Geisler, E. K. Grebel, & D. Minniti, 390 
Hillenbrand, L. A., Strom, S. E., Vrba, F. J., & Keene, J. 1992, 

ApJ, 397, 613 
Hosokawa, T. & Inutsuka, S.-I. 2006, ApJ, 646, 240 
Jackson, J. M., et al. 2006, ApJS, 163, 145 

Lu, F. J., Wang, Q. D., Aschcnbach, B., Durouchoux, P., & Song, 

L. M. 2002, ApJL, 568, 49L 
Malfait, K., Bogaert, E., & Waelkens, C. 1998, A&A, 331, 211 



Martins, F., Schaerer, D., & Hillier, D. J. 2005, A&A, 436, 1049 
MuzeroUe, J. et al. 2004, ApJS, 154, 379 
Nakajima, Y. 2005, AJ, 129, 776 
Onaka, T. et al. 2007, PASJ, 59, S401 

Schaller, G., Schaerer, D., Meyner, G., & Maeder, A. 1992, A&AS, 
96, 269 

Schniidt-Kaler, T. 1982, in Landoldt-Bornstein, New Series Group, 
VL Vol. 2, ed. K. Schaifers, & H. H. Voigt (Berlin: Springer- 

Verlag), 1 

Slane, P. 2007, preprint (arXiv; 0711.1002) 

Spaans, M., Tielens, A.G.G.M., van Dishoeck, E. F., & Bakes, 

E.L.O. 1994, ApJ, 437, 270 
Stil, J. M., Taylor, A. R., Dickey, J. M., Kavars, D. W., Martin, 

P. G., Rothwell, T. A., Boothroyd, A. I., Lockman, F. J., & 

McClure-Griffiths, N. M. 2006, AJ, 132, 1158 
Testi, L., Palla, F., & Natta, A. 1998, A&AS, 133, 81 
Velusamy, T. & Becker, R. H. 1988, AJ, 95, 1162 
Zinnecker, H. & Yorke, H. W. 2007, ARAA, 45, 481 



Star- Forming Infrared Loop around G54. 1+0.3 




Fig. 1. — (left) AKARI 15 fim image of the SNR G54. 1+0.3. The cross marks the position of the pulsar at (0:2000 1 "52000) = 
(19'^30™30M3, +18°52'14."1). VLA 4.85 GHz radio contours are overlaid. The contour levels are equally spaced at every 1.28 K from 
0.22 K in brightness temperature. Contour levels increase toward the center of the SNR. The scale bar corresponds to 1' or 2.3 pc at 8 
kpc. The color bar is given at the bottom of the image in units of MJy sr~^. North is up and east is to the left, (right) Three-color image 
generated from the Spitzer IRAC 5.8 fim (B), the AKARI 15 ^im (G), and the Spitzer MIPS 24 ^m (R). 
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Fig. 2. — Spitzer color-color diagram for 64 MIR sources detected in the IRAC and MIPS 24 fim bands with a [24] magnitude error 
less than 0.2 mag. The sources along the IR loop are well defined in this diagram by their strong [8] — [24] excess and moderate or small 
[3.6] — [5.8] excess, i.e., 5 <[8] — [24]< 9 and <[3.6] — [5.8]< 1.2. The sources within this box are marked by solid dots. (See the text for an 
explanation of half-filled dots.) The areas occupied by class O/I and class II sources determined from previous studies are shown as dotted 
lines. The arrow is a reddening vector of Ay = 10. The inset shows the distribution of the sources within the dashed box on the AKARI 
15 iim image. 
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Fk;. 3. The 2MASS JHKs color-color and color-magnitude diagrams for all 110 MIR som'ces identified in the 2MASS JHKs bands. 
The small dots represent sources not detected in either Spitzer MIPS or IRAC bands while the other symbols have the same meaning as 
in Fig. 2. The dashed line represents the MS. In the color-color diagram, the dotted line represents interstellar reddening line of AO star. 
Along the reddening line, the positon at which Ay = 10 is marked. In the color-magnitude diagram, the dashed line with spectral types 
represents the MS at a distajice of 8 kpc. The dotted lines represent the reddending lines that bound the YSOs. Note that the YSOs (solid 
dots) axe clustered around the positions of reddened B1.5-08 stars. 




Fig. 4. — BED of the most luminous YSO in the Ks band. The open symbols represent observed fluxes while the filled symbols represent 
extinction-corrected fluxes. The Spitzer 70 fim flux is an upper limit. The dashed line is a model SED composed of a stellar emission at 
T* = 35, 000 K and a black body at 190 K. The solid line shows the reddened SED and the dotted lines their individual components. 



